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Abstract 

To support efforts on cooling and trapping of alkaline-earth atoms and designs of atomic clocks, 
we performed ab initio relativistic many-body calculations of electric-dipole transition amplitudes 
between low- lying states of Mg, Ca, and Sr. In particular, we report amplitudes for l P° — > 
1 So, 3 S*i, 1 D2, for ?J P° — > 1 Sq, D2, and for 3 P£ — > 1 £>2 transitions. For Ca, the reduced matrix 
element (4s4p 1 P°||Z)||4s 2 ^o) is in a good agreement with a high-precision experimental value 
deduced from photoassociation spectroscopy [Zinner et al, Phys. Rev. Lett. 85, 2292 (2000) ]. 
An estimated uncertainty of the calculated lifetime of the 353^^° state of Mg is a factor of 
three smaller than that of the most accurate experiment. Calculated binding energies reproduce 
experimental values within 0.1-0.2%. 

PACS numbers: 31.10,+z, 31.15.Ar, 31.15.Md, 32.70.Cs 
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I. INTRODUCTION 



Many-body methods have proven to be a highly accurate tool for determination of atomic 
properties, especially for systems with one valence electron outside a closed-shell core |IJ . For 
alkali-metal atoms a comparison of highly-accurate experimental data with calculations |2| 
allows one to draw a conclusion that modern ab initio methods are capable of predicting 
basic properties of low-lying states with a precision better than 1%. 

For divalent atoms such a comprehensive comparison was previously hindered by a lack 
of high-precision measurements of radiative lifetimes. Despite the lifetimes of the lowest 
nsnp l P° and nsnp 3 P° states were repeatedly obtained both experimentally and theoretically 
H |, |, |, 0, |, 0,00,§00§§ |T7|, Sl^^^^l, persistent 



discrepancies remain. Only very recently, Zinner et al. |3| have achieved 0.4% accuracy for 
the rate of 4s4p l P° — > 4s 2 ^So transition in calcium. This high-precision value was deduced 
from photoassociation spectroscopy of ultracold calcium atoms. One of the purposes of the 
present work is to test the quality of many-body techniques for two- valence electron systems 
by comparing our result with the experimental value from Ref. ||. 

We extend the earlier work |25| and report results of relativistic many-body calcula- 
tion of energy levels and electric-dipole transition amplitudes for Mg, Ca and Sr. The 
calculations are performed in the framework of configuration-interaction approach coupled 
with many-body perturbation theory [^, 27[]. We tabulate electric-dipole amplitudes for 
l P? -> l S , 3 S 1 , 1 D 2 , for 3 Pf -> ^o, 1 ^, and for 3 P 2 ° -> l D 2 transitions and estimate 
theoretical uncertainties. 

Cooling and trapping experiments with alkaline-earth atoms were recently reported for 



Mg H , Ca [||, £9] , and Sr [28|, |3(| • The prospects of achieving Bose-Einstein condensation 
were also discussed || Our accurate transition amplitudes will be helpful in designs 
of cooling schemes and atomic clocks. In addition, these amplitudes will aid in determi- 
nation of long-range atomic interactions, required in calculation of scattering lengths and 
interpretation of cold-collision data. For example, dispersion (van der Waals) coefficient 
Cq characterizes the leading dipole-dipole interaction of two ground-state atoms at large 



internuclear separations . The coefficient C§ is expressed in terms of energy separations 
and electric-dipole matrix elements between the ground and excited atomic states. Approx- 
imately 80% of the total value of C§ arises from the principal transition nsnp l P° — ns 2 1 Sq, 
requiring accurate predictions for the relevant matrix element. Therefore our results will be 
also useful in determination of dispersion coefficients. 



II. METHOD OF CALCULATIONS 

In atomic-structure calculations, correlations are conventionally separated into three 
classes: valence- valence, core-valence, and core-core correlations. A strong repulsion of 
valence electrons has to be treated non-perturbatively, while it is impractical to handle the 
two other classes of correlations with non-perturbative techniques, such as configuration- 
interaction (CI) method. Therefore, it is natural to combine many-body perturbation theory 
(MBPT) with one of the non-perturbative methods. It was suggested p6| to use MBPT to 
construct an effective Hamiltonian H e g defined in the model space of valence electrons. En- 
ergies and wavefunctions of low-lying states are subsequently determined using CI approach, 
i.e. diagonalizing H e g in the valence subspace. Atomic observables are calculated with ef- 
fective operators [0. Following the earlier work, we refer to this method as CI+MBPT 
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formalism. 

In the CI+MBPT approach the energies and wavefunctions are determined from the 
Schrodinger equation 

# eff (£„)|$„> = £ n |$ n >, (1) 
where the effective Hamiltonian is defined as 

H cS (E) = H FC + Z(E). (2) 

Here ifpc is the two-electron Hamiltonian in the frozen core approximation and E is the 
energy-dependent correction, involving core excitations. The operator £ completely accounts 
for the second order of perturbation theory. Determination of the second order corrections 
requires calculation of one- and two-electron diagrams. The one-electron diagrams describe 
an attraction of a valence electron by a (self-) induced core polarization. The two-electron 
diagrams are specific for atoms with several valence electrons and represent an interaction 
of a valence electron with core polarization induced by another valence electron. 

Already at the second order the number of the two-electron diagrams is large and their 
computation is very time-consuming. In the higher orders the calculation of two-electron 
diagrams becomes impractical. Therefore we account for the higher orders of MBPT indi- 



rectly. It was demonstrated |33| that a proper approximation for the effective Hamiltonian 
can substantially improve an agreement between calculated and experimental spectra of 
multielectron atom. One can introduce an energy shift 5 and replace — > H(E — 5) 
in the effective Hamiltonian, Eq. @. The choice £=0 corresponds to the Brillouin-Wigner 
variant of MBPT and the Rayleigh-Schrodinger variant is recovered setting 5 = E n — E^\ 
where E^ is the zero-order energy of level n. The latter is more adequate for multielectron 



systems |34| ; for few-electron systems an intermediate value of 5 is optimal. We have deter- 
mined 5 from a fit of theoretical energy levels to experimental spectrum. Such an optimized 
effective Hamiltonian was used in calculations of transition amplitudes. 

To obtain an effective electric-dipole operator we solved random-phase approximation 
(RPA) equations, thus summing a certain sequence of many-body diagrams to all orders 
of MBPT. The RPA describes a shielding of externally applied field by core electrons. We 
further incorporated one- and two-electron corrections to the RPA to account for a difference 
between the and V^ -2 potentials and for the Pauli exclusion principle. In addition, the 
effective operator included corrections for normalization and structural radiation f2"?fl . The 
RPA equations depend on transition frequency and should be solved independently for each 
transition. However, the frequency dependence was found to be rather weak and we solved 
these equations only at some characteristic frequencies. To monitor a consistency of the 
calculations we employed both length (L) and velocity (V) gauges for the electric-dipole 
operator. 

The computational procedure is similar to calculations of hyperfine structure constants 
and electric-dipole amplitudes for atomic ytterbium |35|, |36[ . We consider Mg, Ca and Sr as 
atoms with two valence electrons above closed cores [ls,...,2p 6 ], [ls,...,3p 6 ], and [ls,...,4p 6 ], 
respectively |37j. One-electron basis set for Mg included ls-13s, 2p-13p, 3d-12d, and 4/- 



11/ orbitals, where the core- and 3,4s, 3,4p, 3,4c?, and 4/ orbitals were Dirac-Hartree-Fock 
(DHF) ones, while all the rest were virtual orbitals. The orbitals ls-3s were constructed by 
solving the DHF equations in approximation, 3p orbitals were obtained in the V^ -1 ap- 
proximation, and 4s, 4p, 3,4d, and 4/ orbitals were constructed in the V^ -2 approximation. 



We determined virtual orbitals using a recurrent procedure, similar to Ref. [38] and employed 



in previous work p6|, [27|, |35], [3(J. The one-electron basis set for Ca included ls-13s, 2p-l3p, 
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3d— 12d, and 4/— 11/ orbitals, where the core- and 4s, 4p, and 3d orbitals are DHF ones, 
while the remaining orbitals are the virtual orbitals. The orbitals ls-4s were constructed by 
solving the DHF equations in the approximation, and 4p and 3c? orbitals were obtained 
in the V^ -1 approximation. Finally, the one-electron basis set for Sr included ls-14s, 2p- 
14p, 3d-13d, and 4/-13/ orbitals, where the core- and 5s, 5p, and 4c? orbitals are DHF ones, 
and all the rest are the virtual orbitals. The orbitals Is— 5s were constructed by solving 
the DHF equations in the approximation, and 5p and 4c? orbitals were obtained in the 
V N ~ l approximation. Configuration-interaction states were formed using these one-particle 
basis sets. It is worth emphasizing that the employed basis sets were sufficiently large to 
obtain numerically converged CI results. A numerical solution of random-phase approxima- 
tion equations required an increase in the number of virtual orbitals. Such extended basis 
sets included ls-ks, 2p-kp, 3d-(k-l)d, 4/-(fc-4)/, and 5g-(k-8)g orbitals, where fc=19, 20,21 
for Mg, Ca, and Sr, respectively. Excitations from all core shells were included in the RPA 
setup. 



III. RESULTS AND DISCUSSION 
A. Energy levels 



In Tables | III] we present calculated energies of low-lying states for Mg, Ca, and 
Sr and compare them with experimental values. The two-electron binding energies were 
obtained both in the framework of conventional configuration-interaction method and using 
the formalism of CI coupled with many-body perturbation theory. Already at the CI stage 
the agreement of the calculated and experimental energies is at the level of 5%. The addition 
of many-body corrections to the Hamiltonian improves the accuracy by approximately an 
order of magnitude. Finally, with an optimal choice of parameter 5 the agreement with 
experimental values improves to 0.1-0.2%. 

Compared to the binding energies, fine-structure splitting of triplet states and singlet- 
triplet energy differences represent a more stringent test of our method. For the 3 P° 2 3 -states 
the fine-structure splitting is reproduced with an accuracy of several per cent in the pure CI 
for all the three atoms, while the 3 P° - 1 P° energy differences are less accurate (especially 
for Ca and Sr). As demonstrated in Ref. |3"3" |, the two-electron exchange Coulomb integral 



Rnp,ns,ns,n P (n=3,4,5 for Mg, Ca, and Sr, respectively) determining the splitting between 3 P° 
and 1 P° states is very sensitive to many-body corrections. Indeed, with these corrections 
included, the agreement with the experimental data improves to 1-2% for all the three atoms. 

The case of the even-parity 3,1 .Dj-states is even more challenging. For Ca, these four states 
are practically degenerate at the CI stage. A repulsion of the level 1 D 2 from the upper-lying 
levels of np 2 configuration pushes it down to the level 3 D 2 and causes their strong mixing. 
As seen from Table || these states are separated only by 10 cm -1 , while the experimental 
energy difference is 1550 cm -1 . As a result, an accurate description of superposition of 
3 D 2 and l D 2 states is important. The 3 D 2 - l D 2 splitting is restored when the many-body 
corrections are included in the effective Hamiltonian. These corrections properly account 
for core polarization screening an interaction between sc? and p 2 configurations. 

For Sr, the fine-structure splittings of 3 Dj states and energy difference between the 3 Dj 
and the 1 D 2 levels are also strongly underestimated in the pure CI method. Again the 
inclusion of the many-body corrections substantially improves the splittings between the in- 
states. It is worth emphasizing, that for such an accurate analysis a number of effects was 
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taken into account, i.e., spin-orbit interaction, configuration interaction, and core- valence 
correlations. A proper account for all these effects is of particular importance for determi- 
nation of electric-dipole amplitudes forbidden in L^-coupling, such as for 3 Pj — > 1 So, 1 D 2 
transitions. 



B. Transition amplitudes 

In this section we present calculations of electric-dipole (El) amplitudes for 3 ' 1 P° — > 1 S , 
3,ipo _^ 1 D 2 , 3 -P 2 ° — y 1 D 2 , and 1 P° — > 3 S\ transitions. The calculated reduced matrix 
elements for Mg, Ca, and Sr are presented in Tables [TV] and |V[ For a transition I —> F the 
Einstein rate coefficients for spontaneous emission (in 1/s) are expressed in terms of these 
reduced matrix elements (i^||£)||/) (a.u.) and wavelengths A (A) as 

2.02613 x 10 18 |(F||D||/)| 2 

Afi = 277T1 - • (3) 

A number of long-range atom-atom interaction coefficients could be directly obtained from 
the calculated matrix elements. At large internuclear separations R an atom in a state \A) 
predominantly interacts with a like atom in a state \B) through a potential V(R) ~ ±C*3/i? 3 , 
provided an electric-dipole transition between the two atomic states \A) and \B) is allowed. 
The coefficient C3 is given by 




C 3 \ = \(A\\D\\B)\* £ (! + W n% „ £L , (4) 



where Q is the conventionally defined sum of projections of total angular momenta on in- 
ternuclear axis. 

From a solution of the eigen- value problem, Eq. ([!]), we obtained wave functions, con- 
structed effective dipole operators, and determined the transition amplitudes. The calcula- 
tions were performed within both traditional configuration-interaction method and CI cou- 
pled with the many-body perturbation theory. The comparison of the CI and the CI+MBPT 
values allows us to estimate an accuracy of our calculations. As it was mentioned above, to 
monitor the consistency of the calculations, we determined the amplitudes using both length 
and velocity gauges for the dipole operator. In general, dipole amplitudes calculated in the 
velocity gauge are more sensitive to many-body corrections; we employ the length form of 
the dipole operator in our final tabulation. 

We start the discussion with the amplitudes for the principal nsnp P° — > ns 2 ^So tran- 
sitions (n = 3 for Mg, n = 4 for Ca, and n = 5 for Sr). Examination of Table [TV] reveals 
that the many-body effects reduce the L-gauge amplitudes by 1.6% for Mg, 5.5% for Ca, 
and 6.4% for Sr. Further, the MBPT corrections bring the length and velocity- form results 
into a closer agreement. For example, for Sr at the CI level the velocity and length forms 
differ by 2.7% and this discrepancy is reduced to 0.8% in the CI+MBPT calculations. 

A dominant theoretical uncertainty of the employed CI+MBPT method is due to im- 
possibility to account for all the orders of many-body perturbation theory. It is worth 
emphasizing that in our CI calculations the basis sets were saturated and the associated 
numerical errors were negligible. We expect that the theoretical uncertainty is proportional 
to the determined many-body correction. In addition, we take into account the proximity 
of the amplitudes obtained in the L- and V-gauges. We estimate the uncertainties for the 
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nsnp l P° — > ns 2 transition amplitudes as 25-30% of the many-body corrections in the 
length gauge. The final values for (nsnp 1 P°\\D\\ns 21 So), recommended from the present 
work, are 4.03(2) for Mg, 4.91(7) for Ca, and 5.28(9) a.u. for Sr. 

We present a comparison of our results for (nsnp l P°\ \D\\ns 2 x Sq) with experimental data 
in Table f[V] and in Fig. [I]. Our estimated accuracy for Mg is a factor of three better than that 
of the most accurate experiment and for Sr is comparable to experimental precision. For 
Ca, the dipole matrix element of the l P° — > x Sq was recently determined with a precision of 
0.2% by Zinner et al. M using photoassociation spectroscopy of ultracold Ca atoms. While 
our result is in harmony with their value, the experimental accuracy is substantially better. 
An updated analysis |40| of photoassociation spectra of Ref. |§ leads to a somewhat better 



agreement with our calculated value. 

A very extensive compilation of earlier theoretical results for the 1 P° — > x 5o transition 
amplitudes can be found in Ref. |J for Mg and in Ref. || for Ca. In a very recent mul- 
ticonfiguration Hartree-Fock (MCHF) calculations for Mg the authors have determined 
(3s3p 1 P°\\D\\3s 21 S ) = 4.008 a.u. This value agrees with our final result of 4.03(2) a.u. 
For heavier Sr the correlation effects are especially pronounced and only a few calculations 
were performed. For example, MCHF calculations for Sr || found in the length gauge 
(5s5p l P°\ \D\ 1 5s 2 ^So) = 5.67 a.u. By contrast to the present work, the core-polarization 
effects were not included in this analysis. As a result, this calculated value is in a good 
agreement with our result 5.63 a.u. obtained at the CI stage, but differs from the final value 
5.28(9) a.u. 

Another nonrelativistically allowed transition is 1 P 1 ° — > 1 D 2 and one could expect that 
this amplitude can be determined with a good accuracy. For Mg this is really so. However, 
for Ca and Sr an admixture of the configuration p 2 brings about large corrections to this 
amplitude, especially in the velocity gauge. Another complication is the following. The 
matrix element of electric-dipole operator can be represented in the V-gauge as (atomic 
units h — \e\ — m e — 1 are used): 

(F\T>\I) = ic(F\ a \I}/(Ej-E F ). (5) 

Here c is the speed of light, Ej and Ep are the energies of initial and final states, and a 
are the Dirac matrices. For the l P° — > X D 2 transition in Ca and Sr the energy denominator 
is approximately 0.01 a.u. Because the El-amplitudes of these transitions ~ 1 a.u. (see 
Table [TV]), the respective numerators are of the order of 0.01 a.u. Correspondingly the matrix 
elements (F\a\I) are small and are very sensitive to corrections, i.e., the V-gauge results are 
unstable. As a result we present only the L-gauge values for 1 P 1 ° — ► X D 2 El amplitudes for 
Ca and Sr. An absence of reliable results in V-gauge hampers an estimate of the accuracy, 
so we rather conservatively take it to be 25%. Note that even with such a large uncertainty 



our value for Sr significantly differs from the experimental value p4| . The measurement in 
|24| has been carried out on the l D 2 — > 1 Sq transition and an interference between electric- 
quadrupole (E2) and Stark-induced dipole amplitudes was observed. In order to determine 
the transition rate a theoretical value of the E2-amplitude for the l D 2 — > 1 S'o transition was 



taken from ETJ . It may be beneficial either to measure directly the rate of the El-transition 



1 P 1 ° — > l D 2 or to measure the rate of the E2-transition 1 D 2 — > 1 S'o. 

For the 3 Pj — > 1 Sq, 1 D 2 transitions the respective El-amplitudes are small; these are non- 
relativistically forbidden intercombination transitions and consequently their amplitudes are 
proportional to spin-orbit interaction. The calculated reduced matrix elements are presented 
in Table |Vf 
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One can see from Tables | ( Q] that the MBPT corrections to the fine structure splittings 
are large, amplifying significance of higher order many-body corrections. In addition, higher 
order corrections in the fine-structure constant a to the Dirac-Coulomb Hamiltonian are 
also important here. As demonstrated in Ref. 0, the Breit interaction reduces the dipole 
amplitude of 3 P° — > 1 Sq transition in Mg by 5%. At the same time for all the intercom- 
bination transitions the agreement between L- and V-gauges is at the level of 6-8%. We 
conservatively estimate the uncertainties of the calculated intercombination El amplitudes 
to be 10-12%. 

To reiterate, we carried out calculations of energies of low-lying levels and electric-dipole 
amplitudes between them for divalent atoms Mg, Ca, and Sr. We employed ab initio rela- 
tivistic configuration interaction method coupled with many-body perturbation theory. The 
calculated removal energies reproduce experimental values within 0.1-0.2%. A special em- 
phasis has been put on accurate determination of electric-dipole amplitudes for principal 
transitions nsnp 1 P° — ► ns 2 1 Sq. For these transitions, we estimated theoretical uncertainty 
to be 0.5% for Mg, 1.4% for Ca, and 1.7% for Sr. For Ca, the reduced matrix element 
(4sAp X P°\ \D\ |4s 2 ^So) is in a good agreement with a high-precision experimental value 
An estimated uncertainty of the calculated lifetime of the lowest l P° state for Mg is a 
factor of three smaller than that of the most accurate experiment. In addition, we evalu- 
ated electric-dipole amplitudes and estimated theoretical uncertainties for 1 P° — > 3 Si, 1 D2, 
3 P° — > ^Sq, l D 2) and for 3 P 2 ° — > X D% transitions. Our results could be useful in designs of 
cooling schemes and atomic clocks, and for accurate description of long-range atom-atom 
interactions needed in interpretation of cold-collision data. 
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TABLE I: Two-electron binding energies E va i in a.u. and energy differences A (cm 1 ) for low-lying 
levels of Mg. 
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Two electron binding energy of the ground state is determined as a sum of the first two 
ionization potentials IP (Mg) and IP (Mg + ), where IP (Mg) = 61669.1 cm -1 and 

IP (Mg+)= 121267.4 cm" 1 @. 



TABLE II: Two-electron binding energies in a.u. and energy differences A in cm 1 for the low-lying 
levels of Ca. 
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a Note that the conventional CI fails to recover the correct ordering of D-states. & For the ground 
state E va i= IP (Ca)+IP (Ca+), where IP (Ca) = 49304.8 cm" 1 and IP (Ca+) = 95752.2 cm" 1 
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TABLE III: Two-electron binding energies in a.u. and energy differences A in cm 1 for the low- 
lying levels of Sr. 
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a For the ground state E va i 
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TABLE IV: Reduced electric-dipole matrix elements for transitions allowed in LS-coupling. n is 
the principal quantum number of the first valence s and p shells and m corresponds to the first 
valence d shell; n = 3 for Mg, 4 for Ca, and 5 for Sr; m = 3 for Mg and Ca, and 4 for Sr. All values 
are in a.u. 
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